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We investigate the effects of the anomalous magnetic moment (AMM) in the equation of state 
(EoS) of a system of charged fermions at finite density in the presence of a magnetic held. In the 
region of strong magnetic fields ( eB > m 2 ) the AMM is found from the one-loop fermion self-energy. 
In contrast to the weak-held AMM found by Schwinger, in the strong magnetic held region the AMM 
depends on the Landau level and decreases with it. The effects of the AMM in the EoS of a dense 
medium are investigated at strong and weak helds using the appropriate AMM expression for each 
case. In contrast with what has been reported in other works, we hnd that the AMM of charged 
fermions makes no signihcant contribution to the EoS at any held value. 
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I. INTRODUCTION 

The fact that strong magnetic fields populate the vast majority of the astrophysical compact 
objects and that they can significantly affect several properties of the star have served as motivation 
for many works focused on the study of the Equation of State (EoS) of magnetized systems of 
fermions and their astrophysical implications [l|-[6|. 

In the presence of a magnetic field B the dispersion relation of charged fermions takes the form E = 
p 2 T 2eBl T m 2 , exhibiting the Landau quantization of the cyclotron frequencies characterized by 
the Landau-level number l = 0,1,2,... [?,]. A magnetic field also affects the density of states which 
now becomes proportional to the field, so the three-momentum integrals change as 



The factor g(l) = [2 — (<Szo)] takes into account the double spin degeneracy of all the Landau 
levels except l = 0. In addition, a magnetic field breaks the rotational SO(3) symmetry, giving 
rise to an anisotropy in the energy-momentum tensor Q and producing a pressure splitting in 
two distinguishable components, one along the field (the longitudinal pressure) and another in the 
perpendicular direction (the transverse pressure). As a consequence, a system of fermions in a 
constant and uniform magnetic field exhibits an anisotropic EoS 0-0. 

Besides modifiying the one-particle Dirac Hamiltonian, a magnetic field can also affect the radiative 
corrections of the fermion self-energy because it introduces an additional tensor F /lv in the theory 
that gives rise to new independent structures like \t(t^ v F iiv , with = §[ 7 ^, 7 ^]- This new term 
corresponds to the coupling between the field and the fermion anomalous magnetic moment (AMM) 
T 0 , which in general can be a function of the magnetic field. It induces a Zeeman splitting in the 
fermion dispersion that removes the spin degeneracy 0 , and the following change in the density of 
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states 


5Zs , ( Z ’ cr )?f^2 [ d P*’ 9(1,0’) = 5io + (1 - Sio) ^ (2) 

l \ n ) J < 7=±1 

with the spin projections a = ±1. 

As shown by Schwinger @ many years ago, at weak fields (eB <C in 2 ), one can make an expansion 
of the fermion self-energy in powers of the magnetic field and find that the leading contribution 
proportional to a^ is linear in the field and given (in natural units) by tB with r = (a/2Tr)nB, 
where hb = e/2m is the Bohr magneton. In this approximation, r is simply independent of the 
magnetic field. At strong field, however, the coefficient of the structure varies as a square 
logarithm of the field 0 and hence cannot be expanded in powers of the magnetic field. At finite 
temperature and/or density, the concept of weak field needs to be revisited, as the field may be 
strong with respect to one of the scales, but weak with respect to the others. 

Charged and massive fermions always possess magnetic moment which in principle can produce in¬ 
teresting physical effects via the modification of the self-energy. Moreover, massless charged fermions 
in the presence of a magnetic field can acquire a dynamical magnetic momen t fl2j -0 through the 
phenomenon of magnetic catalysis of chiral symmetry breaking (MCySB) |17|. The mechanism 
responsible for this effect is related to the dimensional reduction of the infrared dynamics of the par¬ 
ticles in the lowest Landau level (LLL). Such a reduction favors the formation of a chiral condensate 
because there is no energy gap between the infrared fermions in the LLL and the LLL a ntip articles 
in the Dirac sea. This effect has been actively investigated for the last two decades 0-EU- I 11 the 
original studies of the MCySB phenomenon [13]- EH) the catalyzed chiral condensate was assumed 
to generate only a fermion dynamical mass. Recently, however, it has been shown that in QED 
0-0, as well as in quark systems with 0 and without 0 finite density, the MC\SB inevitably 
leads to a dynamical AMM together with a dynamical fermion mass. In massless QED, the dynam¬ 
ical AMM gives rise to a non-perturbative Lande g-factor, a Bohr magneton proportional to the 
inverse of the dynamical mass, and to the realization of a non-perturbative Zeeman effect E3/E3- 

The AMM term in the Hamiltonian changes the energy spectrum of the fermions and can affect in 
principle the properties of the system. Notice that certain neutral particles which, like the neutron, 
are composed by charged particles (charged quarks in this case), can also have nonzero AMM. The 
effects of the nucleons’ and quarks’ AMM on the statistics of magnetized matter have been discussed 
in many works The AMM has been linked among other strong-field effects to stiffening the 

EoS in magnetized stars and to a dramatic variation of the particle fraction, which at very high 
magnetic fields would lead, for example, to pure neutron matter (in the papers of Ref. a a one or 
both of these findings are discussed). However, when investigating the effects of the AMM in any 
physical process we should be careful in using the analytic expression of the AMM that is consistent 
with the magnetic-field strength under consideration. In particular, as we will discuss in detail 
below, considering a linear-in-B approach, which is the approximation used throughout all the Refs, 
a a , is only consistent in the region of weak magnetic fields where a large number of Landau levels 
are occupied because y/eB is smaller than all the other energy scales like mass, temperature, and 
chemical potential. 

The critical field, below which the weak approximation is reliable then depends on the content 
of the stellar matter, the temperature and the density. If one can ignore T and fj ,, the critical 
field separating the weak and strong field regions is determined by the particle mass. For each 
particle species, the critical field can then be obtained by equating the magnetic energy hio c , where 
lo c = qB/mc is the cyclotron frequency in cgs units, to the corresponding rest-energy me 2 . The 
range of critical fields is then quite wide. For example, for electrons, B^ = 4.4 x 10 13 Gauss, for 
quark matter formed by u and d quarks with current masses m u = m<i = 5 MeV/c 2 , it is Bc U,d ‘ > = 
10 2 B^ = 4.4 x 10 15 Gauss, for protons, whose mass is 938 MeV/c 2 , one finds B^ = 1.6 x 10 2 ° 
Gauss, while charged hyperons, which are much heavier, will have a critical field two orders of 



3 


magnitude larger than protons’. At zero temperature and density, a field larger than the critical 
one for that type of particle will constraint them to their LLL. In a system with different types 
of particles, a field may be strong, hence over the critical field, for some of them and weak, below 
the critical, for others, so care must be taken when using weak and strong field approximations to 
consider such subtleties. 

In the present paper we are interested in revising the role of the AMM in the EoS of systems 
of charged fermions, under both weak and strong magnetic fields. This is a due task given that in 
several of the works that have studied the effects of the AMM in the thermodynamical properties one 
can point out several issues in the way the results have been obtained. One of these issues is that the 
strong-field region has been explored inconsistently considering Schwinger’s result for the AMM of all 
the particles, thus ignoring the existence of different critical fields and the fact that the Schwinger’s 
approximation for the AMM breaks down for fields of the order of or larger than the critical one, as 
pointed out many years ago in Ref. E3- Second, when calculating the pressure, some works have 
ignored the existence of a pressure anisotropy 0-0 in a strong magnetic field, so the results were 
obtained basically using a single pressure. Third, some papers neglected the contribution of the 
Maxwell magnetic pressure proportional to B 2 /2 and claimed that the Schwinger AMM produced a 
significant contribution to the statistical quantities, but they concluded this by considering a region 
of fields where the Schwinger approximation not only breaks down, but the magnetic pressure can 
dominate the matter pressure and erase any possible effect of the AMM, as will be shown in this 
paper. Fourth, the pressure of the magnetized vacuum, that is, the contribution that does not depend 
on temperature or density, was also neglected at strong fields where it can be important. To clarify 
all these issues, we shall analyze, through analytical and numerical calculations, the significance of 
the AMM contribution to the main statistical quantities of the magnetized system in the weak and 
strong-field approximations, as well as to the EoS of dense systems with an interest for astrophysical 
applications. With this goal in mind, we shall investigate the weight of each of the participating 
contributions (Maxwell pressure, vacuum pressure, etc.) into the system EoS for different field 
values. As will be demonstrated, our thorough analysis leads us to conclude that, when working 
consistently, the quantum effect of the AMM of charged particles is negligibly small for the EoS of 
the magnetized system at both weak and strong fields. 

The paper is organized as follow. In Sec. [TT] we give the one-loop self-energy of a charged fermion 
system in the presence of a constant and uniform magnetic field using the Ritus’s method [20l |. 
and find the AMM analytical expression for the different LL’s in the strong-field limit. In Sec. 
uni it is calculated the one-loop thermodynamical potential depending on the AMM in the strong- 
field approximation. The result is given as the sum of the renormalized vacuum contribution, 
the contribution at zero temperature and finite density and the thermal contribution. In Sec. 
hyi for the sake of completeness, we calculate the AMM in the weak-field approximation using a 
combination of Ritus eigenfunction method and proper-time representation. In Sec. |Vl we calculate 
the one-loop thermodynamic potential in the weak-field approximation including the AMM found in 
that approximation. The renormalized vacuum contribution and the zero-temperature finite-density 
contribution are presented up to 0((eB) 4 ) order. In Sec. ED we present the numerical results for 
the main thermodynamic quantities, which depend on the AMM, in the weak- and strong-magnetic- 
filed limits. There, we make a thoughtful analysis to determine the significance of the AMM for the 
EoS of strongly and weakly magnetized systems of charged fermions. Finally, in Sec. IVIII we state 
our concluding remarks. We also include four Appendices. In Appendix A, we give details on the 
calculation of the thermodynamical potential in the strong-field approximation at T ^ 0 and /j, ^ 0. 
In Appendix B, we discuss some issues in the calculation of the effective potential at B ^ 0 in the 
Dittrich’s approach. In Appendix C, we derive the Schwinger propagator at B ^ 0, starting from 
the Ritus’s formalism. In Appendix D, the details of the calculation of the thermodynamic potential 
in the weak-field approximation at T ^ 0 and n ^ 0 are given. 
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II. AMM IN THE STRONG-FIELD APPROXIMATION 

The radiative corrections to the magnetic moment of a charged particle in the presence of a 
magnetic field can be found from the one-loop fermion self energy 


S(x, x') = —ie 2 j fl G(x, x')^ v D pv (x - x'), (3) 

G(x,x') denotes the fermion’s propagator in the presence of a uniform and constant magnetic field 
and D flu (x — x') is the photon propagator. 

One can transform the self-energy to momentum space by using Ritus’s approach 

T,(p,p') = J d*xd i yE p (x)T,(x,y)E l pl (y) = (2n) 4 6 ( ' 4) (p - p')H(l)E l (p), (4) 

Index l denotes the Landau-level number; 11(7) = A(sgn(eI3))i5 ,(:) + 1(1 — S 10 ) is a projector that 
separates the LLL (l = 0), with a single spin projection, from the rest (l > 0) with two; 5^ 4 \p — p') = 
S u S(po — Po)S(p 2 — P 2 )$(P 3 ~ p's)- The Ritus eigenfuntions [20| are given by 

Ej,(*) = E ^p(*)A(tr), l'=7°(E^)t 7 0 (5) 

cr=±l 


with 


A(±) 


/ ± *7 1 7 2 
2 


are spin up (+) and down (—) projectors, and 

E+(x) = Aie- l(poa:0+P2a:2+P3x3) A(p), 
E~(x) = Ni-ie- i( ' PoX ° +P2x2+P3x3 ' > Di-i(p) 


( 6 ) 


(7) 


with normalization constant Ni = (AireB ) 1 ^ 4 and Di(p) are the parabolic cylinder functions of 
argument p = \/2eB(xi — p 2 /eB). 

In momentum space the general structure of the self energy is [Ti)| 


s \p) = z\p jf 7 jl + Z l ± p p ± ri + MB + z7I 7 1 7 2 , 


( 8 ) 


Notice the separation between parallel py = (p°,0,0,p 3 ) and perpendicular = (0,0, \/2 eBl, 0) 
components due to the spatial symmetry breaking in a magnetic field that only leaves intact the 
subgroup of rotations along the field direction. In ©, , Z l ± are the wave function’s renormalization 

coefficients. The coefficients Mi and 71 are respectively the radiative corrections to the mass and 
the magnetic moment. Each of them has to be determined as a solution of the Schwinger-Dyson 
(SD) equations of the theory at the given approximation. 

In the one-loop approximation, the Schwinger-Dyson equation leads to an infinite set of couple 
equations that take the form [12b 14] 

E z (p)n(l) = —ie 2 (2eB)Yl(l) [ f Q 6 ^ [L t + L l+1 + L^], 1 = 0,1,2,.... 

J (2tt) 4 q z 


(9) 






5 


with 


and fermion propagator 


L i = llG\p- gbjj, 

Li±i = A (±) 7 ±G l±1 


ip - A(±) 


G\p) = P J 2 1+ ™ W), 


p 2 — m 2 


( 10 ) 

( 11 ) 


Here, we introduced the notation q], = q^/y/2\eB\, p = (p°, 0, — y/2\eB\l,p 3 ) and ( p — q ) M = 

(p° — < 7 °, 0, — \j2\eB\l 1 p 3 — q 3 ). Henceforth, we assume eB >0. As it will become clear below, the 
representation © of the self-energy is particularly convenient for strong-field calculations. 

From Eqs. © and (©, we can extract the equations for the AMM at each LL. In Euclidean space 
they are 


E 0 = (M 0 + To) = e 2 m(4eH) 


d 4 


q e 




(27T) 4 


q 2 


iP ~ o) o + m 2 ip ~ q) i + m 2 


( 12 ) 


/ ~T' 2 

// = — e m 


(2 eB) j ■ 


f d A q e b 

1 

1 

/ (2tt) 4 q 2 ■ 

[ip - q?i+i + m 2 

ip-q) 2 i -1 +m 2 


l > 1 


(13) 


Eq. (fl2l) reflects the single spin orientation of the fermions in the LLL (l = 0) and hence the 
impossibility of determining M 0 and 7o independently mm. Thus, E 0 cannot be interpreted as 
an AMM term, but as the radiative correction to the rest-energy of the LLL particles. Notice, that 
Eq will not produce any Zeeman splitting in the modes of the LLL quasiparticles. 

In the infrared limit po = 0,ps —> 0, and considering the strong-field approximation, the leading 
contributions to ( 1121 ) and (fldl) for l = 1 are respectively given by 


E 0 = M 0 + Tq — 


e 2 m 

8n 3 



e~fi- 1 

q 2 <fjj + in 2 


= m— \u 2 (m 2 / 2eB), 
47T 


(14) 


71 


e 2 m 

16-7T 3 


dq»dq_ 


-q± 


q 2 


% 


= m — ln 2 (m 2 / 2 eH), 
87T 


(15) 


Note that the leading contribution in (fl5l) comes from the spin-down particles in the first LL (the 
term Z — 1 in Cl for l = 1 ), since V2eBl acts, for l > 1 , as a suppressing factor in the denominator 
of the fermion propagator. The result (ThU) coincides with that obtained many years ago in Ref. [I] 
using a different method. As in massless-QED [H|,[l3j], the relation 7i = Eq/2 is satisfied here too. 

For the remaining 71, l > 1, we have 


71= - 


am 

107T 2 


— Mf 


{-7 r [0, 




-e 2 7 F[ 0 ,—M 2 .,] 


— i7re 2 lnM 2 _i — e 2 E i (M/_ 1 ) lnM 2 _i + «7rlnM 2 + i + Ei(M 2 +1 ) lnM 2 + i 


- g; 


LO/ 1,1 _ ]Cf 2 )+P 2 f; 3 ' 0 f 4>1 — m 2 U 

,3^ 0,0,0 lvl l+l) + e Lt 2,3 (^0,0,0 ) / ) 


(16) 


with M 2 ±1 = m 2 + (/±l), 7 ~ 0.577216 is the Euler’s constant, Ei[z\ denoting the exponential integral 
function, r[0,z] is the incomplete gamma function and G////"' (° h \ /'j'/ \ z) the Meijer G-function [2l| . 
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FIG. 1: Comparison between the AMMs Ti and 7i in the strong field region for Landau levels l > l.The plot shows 
a sharp decrease of the ratio 7z>i/7i with increasing Landau levels for two field values. The largest values of 71>i 
occur for the lowest (l > l)-value, but they are still two orders of magnitude smaller than 71- The stronger the field, 
the smaller the values of 71>i for the same Landau level, and the quicker they approach to their asymptotic negligible 
value. 


In Fig.l we show how the AMM’s for l > 1 decrease with respect to 7i as the LL increases. Notice 
that the AMM at strong field is relevant only for the first Landau level, where it grows as the square 
logarithm of the field. As seen in FigllJ already in the second LL the AMM decreases in two orders 
with respect to its value at l = 1 , T 2 /T 1 ~ 0.0668 for rh = 0.1. 

One can explicitly see from (fTHl - dlBl) that, in contrast to the AMM found by Schwinger [oj in 
the weak-field limit, which was the same for all LL’s and had a linear dependence with the field, 
at strong field the AMM does not depend linearly on the field and is different for each LL. Clearly, 
using the Schwinger AMM in the strong-field region would be totally inconsistent and care should 
be taken not to draw any physical conclusions obtained with such a wrong approach. 


III. THERMODYNAMIC POTENTIAL WITH AMM IN THE STRONG-FIELD 

REGION 

To investigate the effects of the AMM in the EoS we consider an effective theory on which the 
fermion propagator is dressed by the one-loop fermion self-energy in the magnetic field, which 
depends on the AMM. The lack of Zeeman splitting (see Eq. (2)) in the LLL separates the propagator 
in the LLL from those in the rest of the levels, so the dressed propagators take the form 

Gq\p) = (P\\ ‘ 7 11 -m)A(+) 


(17) 
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and 


G, 1 (p) = p ■ 7 - m - iTrt 1 7 2 


(18) 


with p = (ip 4 , 0, y/2\eB\l,p 3 ) for l = 0,1,2,... in Euclidean space. Notice that in m we do not 
include the correction Eq, neither in (fl8l) the one-loop corrections to the mass, as they are negligible 
compared to the renormalized mass at B = 0, m. However, since 71 gives rise to a new Lorentz 
structure, it is included in (fl8l) . 

The fermion contribution to the thermodynamic potential of this effective theory is 


niB.p.T, = ~ 


£J 

P* -r. 


TTTTa Indet G 0 1 (p*) 

(2n) z 


E EE j 

<T—±1 1=1 pi 


(^2 hrdet G, V) 


(19) 


Here, /3 = 1/T denotes the inverse temperature, /i the fermion chemical potential and p* = (ip 4 — 
/i, 0, \f2eBl, p 3 ). 

Performing the sum in Matsubara frequencies and calculating the determinants in Eq. m we 
obtain 


n(B,p,T) 



( 20 ) 


which can be rewritten as 

1 a R r°° T? pOO -i 

tt(B,fj,,T) = —-2 / dp 3 E \E V al - lA - 2 “2 / ^E o ln (!+ e' d|Br, ' T!_/i| ) (21) 

/47r pal J ~°° pal P 

In these expressions the sum in the energies E va i include particles/antiparticles (p = ±), up/down 
spin (a = ±), and Landau level l indices. For the LLL, only one spin projection contributes and the 
energy becomes 


Ep,o = ^ 


l = 0, 77 = ±1 


( 22 ) 


For each l ^ 0, the AMM separates the energies of up and down spin (a = ±) as 

E va i = p\]pl + (V 2 eBl + m? + crTi) 2 , l ^ 1, er = ±1, p = ±1 (23) 

Adding and subtracting the vacuum term in (1211) , one can write the thermodynamic potential as 
the sum of vacuum (Q vac ), zero-temperature (fl^), and finite-temperature (Up) contributions, 

Q(B, /x, T) = n vac + + flp (24) 


with 


n vac = n(B, 0 , 0 ) 



(25) 



n^ = n(B,p, 0 ) 


(26) 













0,3 = n{B, 0 ,T) = 
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(27) 


As shown in Section II, in the strong-field approximation all the 77 for l > 1 are very small and 
can be neglected. In this approximation all the energy modes, except for 1 = 1, reduce to the modes 
of the undressed theory. 

In Appendix El we give details on the calculation and carry out the renormalization of the vacuum 
term (12511 in the strong-held region to obtain the following renormalized thermodynamic potential 
in the strong-held approximation 

nl = + n® + of, (28) 

with 



1 

87T 2 



eBs coth(egs) — 1 — 


( eBs ) 2 
3 



In 


2eJ3 

m 2 



(29) 


and 


where 


n?= - eB 


47r 2 

+ e(p-M+) 


j 9(p - to) 

r- 2 - 2 2 , A 4 + \/ A 4 - 2 ™ 2 

/ly /i z — ra z — ram- 

l 

ra 


/ ~2 .2 A 4 + x/m 2 — -^i" 

Vm 2 -M+ -M+ In-L-i- 


In 


2, d+Vd ~M 1 


Mr 


2 


1=2 


r~2 ^72 u 2 i M+ vV ~ | 

Hy/H-Mi -Mi In-—- >, 


^ = - 


27r 2 /3 
eg 
27r 2 /3 
eg 


f dp 3 ln(l + e -/ 3 lv / P§+ m2 +/d)(i + e -/3| Vp|+m 2 -^l) 

Jo 

r dp 3 ln(l + e -/ 3| V / Pi+ M i +2 +Ml)(l + e -/3|V / pI+M+ 2 - / ,|) 

Jo 

On l dps hl(l + e -/3! V^i+Mf 2 +Ml)(l + e -glVpl+A4-r 2 -Ml) 
27T^p Jo 

V [ dps ln(l + e - ^'v /p i +Mi2+ H)(l + 6 -/3 |\/pI+ m ‘ 2 -mI) ; 


eB 


Ml = \J 2eB + m 2 ± 71, M = V^egT+m 2 


(30) 


(31) 


(32) 


The leading vacuum contribution to II in the strong-held approximation (see Appendix [A] for 
details) is then 


aB 2 
67T 

aB 2 

67T 


in ( 

i eg 



^ TO 2 

In ( 

^ eg 


\ 

i TO 2 


47T 2 


47T 2 


, 2 eg „ 

111 5 —h 2 

m 


In 2 

( 2eg V 

2 

‘ 2eg ' 

In —tt + 2 

[v 8 tt / 

r TO 2 / 


TO 2 


( 33 ) 
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where in the second line 71 was evaluated using (fl5l) . and the first term was obtained in [22] 
calculating the effective potential at strong field and neglecting the AMM. 

Here the following comment is in order. The contribution of the AMM to the vacuum part of 
the thermodynamic potential has been previously calculated in [22} with the help of the Green’s 
function. In the strong-field region the result reported in [23| was 



Clearly, there is a discrepancy between (1M1) and (1331) . the origin of which can be traced back to 
some inconsistencies in the treatment followed in |23[ to obtain the strong-field result. On the one 
hand, Ref. @ considered the AMM found by Schwinger and used it in calculations at arbitrary 
field strength, including the strong field region, despite that as already discussed, in this region the 
AMM becomes a very different function of the magnetic field and depends on the LL. On the other 
hand, as described in Appendix B, several steps of the calculations done in [22} were only valid for 
weak fields, however, they were used indistinctly for weak and strong fields. The calculations of Ref. 
[23} are then reliable in the weak-field region, but fail to describe the AMM-dependent terms in the 
strong-field case. 


IV. AMM IN THE WEAK-FIELD APPROXIMATION 


To get the AMM in the weak-field limit, it will be convenient, for the sake of clarity and to shade 
light in our discussion, to use an alternative method that combines the Ritus’s approach, where 
the LL contributions are explicit, and the proper time formalism. The result has to coincide with 
that found by Schwinger |2J] by using an independent method. In doing this, we will stress the 
steps where the weak-field approximation becomes a basic element of the derivation. Also, it will be 
apparent the different behavior of the AMM’s for the different LL’s in each approximation. 

Due to the fact that the energy separation between consecutive LL’s is proportional to \/2eB , 
it is expected that at weak field all Landau levels will contribute on equal footing to the fermion 
self-energy. Then, we will find convenient, in order to take into account the contributions of all the 
LL’s into the fermion self-energy, to work with the self-energy operator in the configuration space 
© with the field dependent fermion Green’s function given in the Ritus’s approach as 

G{x,y) = ^(^ 4 Mx)G l m P (y )• (35) 

where 

G\p) = -L- H(0 (36) 

p — m 

is the fermion propagator in momentum space and we introduced the notation 




dp 0 dp 2 dp 3 

(2nr 


We can rewrite the electron propagator in (1551) as 


(37) 


G{x,y) 


(P* + TO )^ 


d A p E p (x)U(l)E p (y) 

(27r) 4 p 2 — m 2 


(38) 


where n x M = idj‘ — eA M and we used the property, E p (x) = E P (x)p, satisfied by the Ritus 
eigenfunctions, E p (a:), defined in ©-©. 
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Now, to perform the summation over all Landau levels, we use the proper-time representation 
and the integral representation for the parabolic cylinder functions, so, we get (See Appendix [Cl for 
details) 


G(x,y) = -(fi x +m) 


(47r) 2 


x e 


ieBs 


A(+) 


JO 

—ieB, 


° eBds 
s sin (eBs) 

'A(-)] 


0 —is(m 2 -ie) -»[^ (x-y) j)- 4 ta °g Ba) 0~3/)l] 


(39) 


where 4>(x, y) = exp \i^§-{x 2 — J/ 2 )(iCi + Vi)\ is the well known Schwinger’s phase (recall that F 2 1 = 
B) H. 

In the weak-field approximation (eB -C m 2 ), we can perform a Taylor expansion up to linear 
terms in eB in the integrand of Eq. (1591) , to find 


G(x,y) 


(4tt) 2 Jo * 2 

x [1 + ieBs( A(+) — A(—))], 





vY 


(40) 


where we used the fact that the Schwinger phase satisfies the identity 

n X li <f>{x,y) = P -F^{x~ y) v $(x,y). 


(41) 


With the help of the identities eB( A(+) — A(—)) = — #< 7 ^F^ and [y p , a^] F pv X p = 
—4i r y pj Fp P X p with X p an arbitrary four-vector, we rewrite Eq. (14CT1) as 


G(x,y ) 


H^,y) r ds is(m 2_ ie) 

(4tt) 2 Jo * 2 



1 

4 




1—0 r 
2 


' F, 


His 



1 


si^cr^F^ 


_£0ws/) [°° ds_ is{ m 2 -ie) p -i<^- 

(47r) 2 Jo s 2 

x _ i ^ + m — msi^a pl 'FpJ^j 


(42) 


where in the last line we kept up to linear terms in eB, in agreement with the weak-field approxi¬ 
mation. 

To find the fermion self-energy we have to substitute in Q the photon propagator in configuration 
space 


Dpvix - y) 


dt 

(47T) 2 Jo t 2 6 


together with the fermion propagator (1551) . Then, after using the identities 7 P 7 P 7 M 
7 At 7 a 7 /3 7 /i = Ag a/3 and 7 p 7 p 7 a 7 /3 7 P = — 2'yP'y a 'y p , the fermion self-energy reduces to 


£ (x,y) 


A e 2 $(x,y) dsdt_ is{m 2 _ ie) i{1+ i 

(4tt) 4 J 0 S 2 t 2 


Now, we introduce a new variable defined as 


(43) 

- 2 7 p , 


(44) 


1 

w 


1 

- + 

s 


1 

7 


(45) 
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to rewrite Eq. (PHI) as follows 


T,(x,y) = i 


. e 2 ${x,y) [°° ds f s dw is(m 2 _ ie) 

— >1 - | - I p v ' g 4 it; 


(4tt) 4 Jo s 2 


w 4 


x - 


t-i i 


2 {^- ^ +4m 


On the other hand, from Eq. m we have 


1 


P x - ^ [#- # , i^cr^F ^J ) 4>(x,y)e l< ^ = 4>(a’, y)e l< 




2w 


(46) 


(47) 


Thus, solving this self-consistent equation iteratively in the weak-field limit, the leading terms 
have the form 


, o»-«r 




d>(x,y)e *' [^f — jf| ~ 2w ( fl x — — 2wfl x , i—a^F^ ) <fr(x,y)e . (48) 


Replacing Eq. (H51) into Eq. flU . we obtain 


T,(x,y) ~ i 


ds 


(47r) 2 
1 

"2 


dw e -i«(™ 2 -^) ( _ 


2w 


(r 


w 

"2 L 


r, *2^^, 


{2wp x , i^F^+Am^j < a:|C/|2/ >, 


(49) 


where we used that the space-representation of the proper-time evolution operator U = e lHw , with 
U = (n M ) 2 , is given by 


< x\U\y >= 1 $(x,y)e l ~ 4 ^ e . 

(47r) z w; z 


(50) 


Multiplying Eq. (l49l) by ip(y), a solution of the Dirac equation, [p — = 0, and integrating 

over y, we arrive at 

„2 /-oo Jo r s / 2w 


d yS(x, y)J>(y) = i 


(4t r) 2 


—w 


{r, 


50 ds 


/o 


dw e" is(m " ie) ( - 
| -(- 4mj 


s 

iw^CTauF^ 


(r - f [r 

J)(x)e iwm2 


i\^ v F, v 


where we used the fact that 

J d A y < x\U\y > J>(y) = J dJy < x\U\y X y\ijj >=< x\U\J> >=< x\J) > e lwm 
Bearing in mind the identity [A, B] = {4, B} — 2BA, then Eq. (I5T1) becomes 

J d 4 yV(x,y)i>(y) = ^ e -^-u) (4m - ^Jl x + tmwi^F^ 

~ w (l + j) {r,i^F^}) p(x)e™ m2 

Next, using once again that fhj = rmp in (l53l) . we have 


(51) 

(52) 


(53) 


d 4 y^(x,y)i/j{y) = i 


dsdw 


(4tt) 2 


j—is(m 2 —ie) 


2m ^2 —+ 2 mw ^1-- j i—a^F 111 ' ip(x) 


(54) 
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where we kept only the leading contribution in the magnetic field O(eB). 
Eq. (l54l) can be written as 


T,(x)ip(x) = (ml — Ti r ) lr ) 2 )ip(x) 


(55) 


Here, the two terms corresponding to the two independent Dirac structures, are the radiative mass 
and AMM, respectively given by 


to = too 


T = 


TO 0 + “ ^ f dw e -i(s- W )^-ie) 

2tt Jo s 2 J 0 \ sJ 

„ pC 

2tt Jo 


dsdw 


0 —i(s—w){rn 2 —ie) 


W ( W\ _ 

— (1-) eBrri 

s V s J 


(56) 

(57) 

(58) 


with too the electron bare mass. These two equations coincide with Schwinger’s results for the mass 
and AMM Ref. 0L 

Following Ref. [241 ] . the last two equations can be evaluated by using a new variable u = 1 — w/s 
and making the replacement s —> —is, so, we arrive at 


TO 

r 


3 a 

mo + -—to 

4 7 T 

a eR 
27 t 2to 



(59) 

(60) 


where to is the electron’s renormalized mass. 

As we can explicitly see in the derivations of this section, the formula (1601) for the AMM is only 
valid in the weak-field limit (eB < to 2 ). 


V. THERMODYNAMIC POTENTIAL WITH AMM IN THE WEAK-FIELD LIMIT 


In this Section we aim to find the thermodynamic potential of the effective theory in a weak 
magnetic field. Since the field is weak, all the LLs contribute to the AMM, so the AMM entering in 
the effective theory is now the one found by Schwinger many years ago and given by the equation 
(1601) . It is well known that the Euler-Heisenberg Lagrangian (i.e. the vacuum contribution to 
the thermodynamic potential) for charged fermions in the weak-field limit contains divergences 
associated with the renormalization of the field and charge [24j • A peculiarity of the present case is 
the appearance of new divergencies connected to the AMM. 

Following the same approach of Section III, we can write the renormalized thermodynamic po¬ 
tential in the weak-field limit as the sum of the vacuum (f2™£^), zero-temperature (flj^) and 
finite-temperature (fig ') contributions 


fi^ = n w ^ + fi w + Q^, 

R vac ' u u u p 7 


(61) 


where 


^ R) = “2 jky /°° J e - 5 K+ r2) |-|eR| S sinh( 2 S rTO) 


+ 


2 ( 2 - 7 t ) 2 ° 3 

ypK n +i °° dt t+ ^r 


2?n J'y — ioo t 2 


e t+ 4 |eR|scoth ( s 1 — 


sT 2 

t 


eB 


- [ i + s r 2 + 2 S 2 TO 2 r 2 -2TO S 2 reR-i s 2 r 4 + ^^ 

6 


3 


( 62 ) 
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Q ^ = - (2^ |(/^-e++o)0(M-e++o) + X! 5 Z(m-£+<t i)6 , (M-e +( ri)| , (63) 

and 

J dp 3 ln(l+e-M e ++°+^ (l+e~M e ++°-^ 

- / dp 3 ln(l + e _/3|E +^ +Ml ) (l + e -/3|£+^-Mlj ? ( 64 ) 

' 7I '^ ^ (T=±l i = l ^ 

with 

Eijcri = ^ + (\/m 2 + 2lcB + crT^ (65) 

and T given in (l60l) . The details of these calculations are given in Appendix D. 

Note that in the above equations the AMM does not contain any dependence on the LL’s. This 
is one of the main differences with respect to the strong-field case. 

In the weak-field expansion, the leading terms of the renormalized potential (1621) . taken up to 
0{{eB ) A ) order and a 2 correction, are (see Appendix D) 


nZi R) (B, o,o) 


1 f(e£) 4 

2(2tt) 2 \ 45m 4 


(eBfT 2 \ 

3m 2 j 


( 66 ) 


Here, we should mention that if in the effective potential found in Ref. [23| we take the weak-field 
limit up to the same order we are considering here, we obtain the same result, although in the 
calculation carried out in Ref. [Isj there are certain inconsistencies as we discuss in Appendix B. 
This is indicating that the inconsistencies in the calculations of 0] did not affect the results in the 
weak-field limit, but only those at strong field. 

Similarly, the leading terms in the weak-field approximation up to 0((eB) 2 ) order and a 2 correc¬ 
tion of the finite-density, zero-temperature thermodynamic potential reads 


= ~g~2 | i^ p3p ~ m2flpF + w4 ln ( PF ?? t ~ ^ 

-S («"' + "* 2 in { EL ^ t ) ) ^ - m)+ in s<fi - m > 

(67) 


where pf = \fpF-m 2 . In the above equation the first three terms are the contribution to the 
thermodynamic potential at B = T = 0 of a system of fermions at finite density j2(| • The next three 
terms are the AMM contributions to the magnetized fermion thermodynamic potential at finite 
density and zero temperature. Finally, the last term accounts for the pure magnetic contribution. 
Comparing (RICH) and (1U71) . we can see that at weak field (eB <C p 2 ), the vacuum contribution can be 
neglected with respect to the finite-density contribution at zero temperature. 


VI. EQUATION OF STATE WITH AMM 

In this section we will investigate the effect of the AMM in the EoS. We will separate the analysis 
for the weak-field and strong-field regimes using the expressions of the thermodynamic potential at 
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T = 0 found in each limit, as well as the corresponding formula of the AMM. Also, to show the 
consequences of using the wrong prescription of taking the AMM in the Schwinger approximation 
when one is dealing with strong-field calculations, we will find the deviation in the parallel pressure 
produced by using T in place of 7i. 


4 

10 



FIG. 2: (Color online) Vacuum effects on the parallel and transverse pressures. Profile of the parallel and transverse 
pressures versus B /in the strong-field region, for p = 10.0 MeV and m/p ~ 0.051. All the plots include the AMM, 
but ignore the Maxwell pressure. Plots with vacuum contribution included (black) and ignored (red) are compared. 
The figure shows that the vacuum terms in the two pressures become progressively relevant with increasing magnetic 
field. In the region of very strong fields the vacuum terms are big enough to produce a sizable negative contribution 
to the perpendicular pressure. 


As pointed out in the Introduction, the breaking of the rotational symmetry produced by a 
uniform magnetic field gives rise to an anisotropy in the energy-momentum tensor j§} that leads 
to a splitting of the pressure into two different components 0]-@i one along the field (longitudinal 
pressure, Fy) and another in the perpendicular direction (transverse pressure, P±). Those pressures 
and the energy density are respectively given by 

P ll= -nf-B 2 / 2, P ± =-nf-BM + B 2 /2, e = fif +/jJ\f + B 2 / 2 (68) 

Here, fij? represents the fermion renormalized thermodynamic potential, given below in the strong- 
field, fig , and weak-field, fi^, approximations; Ai = — (dfl^/dB) is the system magnetization; J\f = 
— (dfi^/dfj.) its particle number density; and the terms proportional to B 2 /2 are the renormalized 
Maxwell contributions @, [271]. In the following considerations, we consider the zero-temperature 
limit, which is of particular interest for astrophysical applications. 
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FIG. 3: (Color online) AMM effects on the parallel and transverse pressures. Profile of the parallel and transverse 
pressures with the field in the strong-field region including vacuum contribution, but ignoring the Maxwell pressure. 
The values p = 10.0 MeV and m/p ~ 0.051 were used. Black (red) curve indicates AMM yt 0 (AMM = 0). The 
zooms in the right panels show that the AMM makes a negligibly small contribution to each pressure in the strong-field 
regime. 


A. EoS in the strong-field approximation 


The zero-temperature limit of the thermodynamic potential in the strong-field limit, fig, is de¬ 
termined by the first two terms of (l28l) . which are given by Eqs. (l33l) and (1301) . respectively 


fl s W(T = 0) = {i%W + n0 = -¥?£ In 


24tt 2 


eB 

rr.2 


^ rj-2 

1 " 


eB 

47T 2 


j^O- 

m) 

/ 2-2 2 , M + vV m 2 

fiy fi z — m z — m in- 

m 




Jy? - M+ 

-M+) 


^ u 2 \f +2 M +2 In ^ + ^ 


M+ 





-Mr) 


!..2 \f~ 2 M~ 2 111 ^ ^ 

/ M 2 - Mi 

fJj IVJ IVJ 111 

M// 


(69) 


To ensure that only the zero and first LL’s are populated for the considered magnetic-field values, 
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FIG. 4: (Color online) AMM effects on the energy density. Energy density in the strong-field region including the 
vacuum contribution and ignoring the Maxwell one, for fi = 10.0 MeV and m/fi ~ 0.051. Black curves indicate 
AMM 7^ 0, red curves indicate AMM = 0. The zoom shows that the AMM makes no significant contribution to the 
energy density. 


we should satisfy the condition m < /i < \/4 eB + m 2 . Using the thermodynamic potential (1691) . 
we find from (1681) the parallel and transverse pressures in the strong-held approximation. Their 
graphical representations as a function of the magnetic held are given in Figj2j without the Maxwell 
term. There, we compare the cases where the vacuum term is included and ignored. As can be seen, 
the vacuum term in the strong-held regime produces a remarkable contribution leading to a quicker 
turn of the transverse pressure toward negative values. 

In Fig. Owe compare the parallel and transverse pressures, with and without the inclusion of the 
AMM. In both cases the vacuum and medium contributions are included, but the Maxwell term 
is ignored. As can be seen in the zooms of the right panel, the AMM in the strong-held regime 
makes a negligibly small contribution to the pressures, since it enters as an a-correction in the 
thermodynamic potential. 

The transverse pressure in Figs J5]and [3] displays a de Haas-van Alphen oscillation associated with 
the transition of the fermion filling from the Landau level l = 1 to the LLL as the magnetic held 
increases. De Haas-van Alphen oscillations are also present in the energy density (Fig.Q) and in 
the magnetization (Figj5|). The size of the de Haas-van Alphen oscillations are not large enough in 
the parallel pressure to be observable at the scales under consideration. It is also evident from FigsQ] 
and [5] that the AMM makes no signiheant contribution to the energy density or the magnetization. 

The effects of the Maxwell term on the pressures in the strong held region are explored in FiglHJ 
Without the B 2 /2 term, the parallel pressure increases with the held (see Figs|3]and[3]). However, 
once this contribution is added, it wins over the rest of the terms in the thermodynamic potential 
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FIG. 5: (Color online) AMM effects on the magnetization. Magnetization versus magnetic field in the strong-field 
region including the vacuum contribution and ignoring the Maxwell one, for fi = 10.0 MeV and m/fi ~ 0.051. Black 
curves indicate AMM ^ 0, red curves indicate AMM = 0. The zoom shows that the AMM makes no significant 
contribution to the magnetization. 


and the field-dependence is inverted, eventually becoming negative for strong enough fields. For 
the transverse pressure, the opposite effect occurs. On the other hand, the Maxwell contribution 
completely erases any difference produced by the AMM among the parallel and transverse pressures. 
In a similar fashion, a Maxwell term modifies the profile of the energy density with the field in the 
strong-field region, making it to increase with the field as seen in Fig0 All these results underline 
that once the field becomes of the order or larger than the chemical potential, the Maxwell term 
dominates the behavior of all the thermodynamical quantities. 

Now, to show the implications of erroneously using the Schwinger formula for the AMM ((501) in 
the strong-field approximation, let us substitute 71 in (l69l) with T = {a/2n){eB/2m) and find the 
parallel pressure including the vacuum term, which we denote as P\\{Schw). Let us define now the 
ratio 


A sch W = \P\\(Schw) — P\\(AMM = 0)| 

A amm \P ll (AMM)^P ll (AMM = 0)\ [ j 

where Fy {AMM) is the parallel pressure defined from (l69l) with 71 taken from (fl5l) . and Fy {AMM = 
0) with 71 = 0. 

From Fig[8]we see that the net effect of erroneously using the Schwinger AMM in the strong-field 
approximation is to change in one order of magnitude the contribution of the AMM to the parallel 
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FIG. 6: (Color online) Maxwell term effects on the pressures. Profile of the parallel and transverse pressures versus 
the magnetic field in the strong-field region including vacuum and Maxwell contributions and using p = 10.0 MeV and 
m/fi ~ 0.051. Black curves indicate AMM ^ 0, red curves indicate AMM = 0. The Maxwell term totally changes the 
profile of the pressures, as can be seen by comparing these plots with those shown in Figs. © and ©. At sufficiently 
high B, the Maxwell contribution, which enters with different signs on each pressure, has a magnitude large enough 
to produce a noticeable decrease (increase) in the parallel (transverse) pressure. The zooms in the right-hand panels 
show how the Maxwell pressure practically erases the already small differences between the pressure with and without 
the AMM. 


pressure. That is, although the Schwinger’s AMM makes an a contribution into the thermody¬ 
namic potential, its linear dependence on the field produces a larger effect at strong field than the 
square logarithmic behavior found for the AMM in the strong-field approximation in (1151) . The kink 
appearing in Fig(8] corresponds to the field value where the de Haas-van Alphen oscillation takes 
place in FigsJS] and [3] Thus, for the scale used in Fig. [8] the de Haas-van Alphen oscillation of the 
parallel pressure becomes apparent. The effect of the Schwinger’s AMM is to stiffer both pressures. 
Nevertheless, when the Maxwell contributions are added one can see that at such high fields the 
effect of the Schwinger’s AMM also becomes negligible. In this analysis we were considering only the 
parallel pressure to make contact with other results in the literature, because the parallel pressure 
is what has been usually considered when the pressure splitting has been ignored, but if we repeat 
the calculation for the transverse pressure the result will be identical. 









19 


10 6 



FIG. 7: (Color online) Maxwell term effects on the energy density. Energy density versus the magnetic field in the 
strong-field region including vacuum and Maxwell contributions and using p, = 10.0 MeV and m/fi ~ 0.051. Black 
(red) curve indicates AMM ^ 0 (AMM = 0). The Maxwell term completely changes the behavior of the energy 
density with the field, as can be gathered by comparing these plots with those of Fig. q. The zoom shows that there 
are no substantial differences when the AMM is present. 


B. EoS in the weak-limit approximation 


In the weak-field approximation, the thermodynamic potential at zero temperature is given by 
the sum of the renormalized leading vacuum term (1661) and the medium contribution (1671) 


QW(R) (T _ o) _ n W(R) , - — I ( eB ^ 

1 Uj + ^“8^145m 4 + 3m 2 J 


-I- 

8tt 2 | 3 


pp 3 F — m 2 fipF + m 4 In 


Pf + V 
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r 2 

47T 

\eB\mT 


( WF +ln 


Pf + M 


TO 


27T 2 


ln 


Pf + P 


TO 


>0{p — to) 


9(p — to) 


or ^ \ eB \ 2 ^ fPF+P 


9{p — to) (71) 


Here the AMM T is taken in the Schwinger approximation (l60l) . 

We can now use m to calculate the parallel and transverse pressures in the weak-field region. We 
plot them, ignoring the Maxwell term, as functions of the magnetic field in Fig|9] For comparison, we 
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FIG. 8: Ratio of the variation of the parallel pressures with AMM’s a la Schwinger and in the strong-field limit, 
A Schw/^AMM — |P|| (Schw) — P|| (AMM = 0)\/\P ll (AMM)-P ll (AMM = 0)|, for magnetic fields in the strong-field 
regime. Notice that the consequence of erroneously using the Schwinger AMM in the strong-field approximation is to 
change the contribution of the AMM to the parallel pressure in one order of magnitude. The kink corresponds to the 
field value where the transition of the fermion filling from the first LL to the LLL takes place. 


also plotted the corresponding pressures with no AMM. It can be gathered from the figure that the 
AMM fails to produce a significant effect on the pressures also in the weak field case (B < B*/ 10). 
This is easy to understand since at weak field, T ~ a(eB/2m) makes a small correction to the 
thermodynamic potential. In contrast to the strong-field regime results, in the weak-field range 
(i eB < m 2 < /x 2 ) the two pressures remain positive. 

Similarly to what happens at strong fields, once the Maxwell pressure is included any effect of the 
AMM is erased and the dependence of the pressures on the magnetic field is inverted with respect 
to the case without it. As shown in Fig [TUI with the Maxwell contribution inserted, the transverse 
pressure increases with the field, while the parallel pressure decreases. However, the latter decreases 
slowly and does not become negative. 

We can express the percentage splitting between the parallel and transverse pressures in the 
weak-field region as 


A[%] 


P± ~ P\\ 1 

\P\\(eB~0) | X 


(72) 


All the quantities in this formula contain the Schwinger AMM denoted by T. As can be gathered 
from FigslU and El at very weak fields (B <C B*) the two pressures coincide, so we have denoted 
them in this region by a common symbol P\\{eB ~ 0)). 
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FIG. 9: (Color online) AMM effects on the parallel and transverse pressures in the weak field region without the 
Maxwell term. Profile of the parallel and transverse pressures with the magnetic field in the weak-field region including 
vacuum contribution, but ignoring the Maxwell pressure. The values p = 10.0 MeV and m/fi ~ 0.051 were used. 
Black (red) curve indicates AMM ^ 0 (AMM = 0). Notice that the AMM does not produce any significant effect 
on the pressures. In the region of interest the two pressures remain positive. 


A plot of the percentage splitting (1771) as a function of the magnetic field is given in FigfTTl using 
= 10.0 MeV. Here we included the Maxwell contribution. We can see that for fields up to B ~ 
the splitting is much smaller than 2%. For strong fields it can be ten times larger [6j. Thus, we 
conclude that in the weak-field regime the pressure splitting is negligibly small and the inclusion of 
the AMM makes no difference. 


VII. CONCLUSIONS 

In this paper we presented a thorough study of the EoS of charged fermions endowed with AMM 
in a dense medium. The calculations were done in the strong and weak field limits, including on 
each case, the corresponding AMM previously obtained in those approximations. 

To obtain the AMM in the strong and weak field limits we took advantage of the Ritus’s approach, 
which allows diagonalizing the self-energy in momentum space and separating the different LL 
contributions. We found that at strong field, the particle gets an AMM that depends on the LL’s, 
with a noticeable difference between the value of the AMM for l = 1, which is proportional to the 
square logarithm of the magnetic field and equal to half the rest energy of the LLL, and the AMM 
for the remaining levels l > 1, which is much smaller and decreases very rapidly with the field. 
In sharp contrast, in the weak-field approximation the AMM is independent of the LL and grows 
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FIG. 10: (Color online) Effects of the Maxwell term on the parallel and transverse pressures in the weak-field region. 
Profile of the parallel and transverse pressures with the field in the weak-field region including vacuum and Maxwell 
contributions. The values p, = 10.0 MeV and rnf // ~ 0.051 were used. Black (red) curve indicates AMM f 0 
(AMM = 0). Due to the Maxwell contribution the transverse pressure increases with the field, while the parallel 
pressure decreases, but it does not cross zero. As in the strong-field case, the AMM effects are totally erased by the 
Maxwell contribution. 


linearly with the field. 

We calculated the different contributions to the thermodynamic potential: vacuum, medium and 
Maxwell terms, and found how they affected the energy density and the parallel and transverse 
pressures in the presence of a magnetic held. We then investigated the role on the EoS of each of 
these contributions, with and without the inclusion of the particle’s AMM in the strong, as well as 
in the weak-field approximations. 

Our main conclusion is that the AMM of charged fermions makes no difference in the EoS, neither 
at strong, nor at weak magnetic fields. This statement contradicts some claims in the literature 
i about the significant effects of the AMM at strong fields. We believe that the origin of the 
discrepancy is connected to using Schwinger’s result to describe the AMM of all the fermion species 
entering in the theory and assuming this is a good approach for arbitrary values of the magnetic 
field. Such an inconsistent treatment to study the effects of the AMM on the thermodynamics of 
the system is a common element of several works of Ref. [3], so the claims about the important 
role of the AMM in the EoS should be taken with a grain of salt, at least for the cases where the 
claims are connected to the role of the AMM of charged fermions, since our results, and hence our 
criticism, applies only to the effects of the AMM of charged fermions. It would be interesting in the 
near future to undertake a similar study about the effect of the AMM of neutral composite fermions. 
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FIG. 11: Profile of the pressure splitting percentage versus the magnetic field. Pressure splitting (Eq. (1721) 1 vs 
magnetic field including the Maxwell pressure for // = 10.0 MeV. Notice that the splitting percentage is smaller than 
2 , so in the weak-field regime the pressure splitting is negligibly small. 
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Appendix A: Thermodynamic potential in the strong-field approximation 

1. Vacuum contribution 

We are interested here in the vacuum contribution to the thermodynamic potential in the strong 
field region fl® ac . It is obtained from il vac in (El) , after taking into account that at strong field all 
7/ with l > 1 can be neglected and hence the energy modes E v<7 i become 


E ±+o = ±\jp\ + m 2 , E ±al = ±\Jpj + (\/m 2 + 2 eB + oTi) 2 , E ±a{l>1) = ±\Jp\ + m 2 + 2 eBl 

(Al) 
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Therefore, 
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1=2 


If we add and subtract 2y / p| + to 2 + 2ei3, it can be written as 
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Performing the sum in LL’s, integrating in momentum, and introducing the proper-time repre¬ 
sentation, we obtain 
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The terms without the coth(eBs) can be integrated using the formula 


— ir-dx = - 1 - aEi{ax ), 


(A4) 


(AS) 


where Ei(ax) is the exponential integral function. Evaluating in the limits of integration and using 
the formula 


Ei(ax) = In \ax\ + 
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which is valid for \ax\ < 1, we obtain that the leading contribution in the limit of large A is 
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The term proportional to coth(eI3s) in (1A4I) has ultraviolet divergencies that can be identified by 
expanding it in powers of s. They reduce to two terms, 


1 

87T 2 


ds 


and 


T/a 2 


ds 


8tt 2 J 1/a 2 s 3 


: ( eBs ) 2 
3 ’ 


(AS) 


(A9) 
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To eliminate these divergences, one can add and subtract (IASI) and (IA9I) to (IA4I) p3 |. The divergence 
(IA8D . which depends on the renormalized mass, is then incorporated into a renormalized vacuum 
constant, while m is absorbed into the bare Maxwell energy to define renormalized field and 
charge. Then mass, charge and magnetic field are all renormalized in the standard way. The second 
term in the rhs of (1A71) is divergent though. Keeping in mind that the parameter 71 is finite, since it 
is given by the AMM found from the renormalized one-loop self-energy of the original theory, hence 
was introduced in the effective theory as a given finite parameter, we can treat this term in the same 
way done in [24] with (IA8I) and incorporate it into the renormalized vacuum constant. 

Therefore, the renormalized vacuum contribution to the thermodynamic potential in the strong 
field limit takes the form 


n S(R) = 


A f 

87r Jl/A 2 


- [eBs coth(eBs) — 1 — 


( eBs ) 2 


eB 

47T 2 


T? 


\ 2 eB 
In —— + 2 


(A10) 


To extract the strong-field contribution to the integral in (IA10I) . one separates it in two pieces 



(All) 


The first integral can be done expanding the coth in powers of its argument and gives a negligible 
contribution in the limit A —> oo. The second integral is finite. Its leading contribution is given 
by In (^), which coincides with the strong-field result found in [29]. Then, keeping just the 

leading-in-a terms in the strong-field approximation, the renormalized vacuum potential reduces to 


n S(R) _ 
vac ~ 67T 


eB q -2 

47T 2 ' 1 


In 


2 eB 


(A12) 


Substituting 7i with the result (l28l) . we obtain 


n S(R) = _^l ln (2l 


6n 


m* 


eB 

47T 2 



In 


2 eB 
m 2 



(A13) 


so clearly the leading vacuum contribution is mainly determined by the strong-field contribution at 
zero AMM [22|], given by the first term in the rhs of (IA13I) . 


2. Finite-density contribution 

At zero temperature the contribution of the fermions in the Fermi sphere to the thermodynamic 
potential is given by 



(A14) 


with energy spectrum given in (IA1I) . 
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Summing in 77 we get 

eB r°° 

^ = ~ J dp 3 [\E +a i — n\ — 2E +a i + ( E +<t i + fi )] 

<J l °° 

^ jj _^ noo 

= - g ^2 ^ / d P3 i E +<rl ~ V + M ~ 2£ +cr z] 0(B +(T z - /i) 

<rl 00 

/ oo 

dp 3 [A 1 — E +cr i + E +a i + p — 2E +cr i\ 0(p — E +a {) 

at -°° 

/ OO 

dp 3 (77 - B +CT ;) 0 (// - B +CT z) 

-00 


eB 

87T 2 


eB 

4-7T 2 


(A15) 


(T l J -°° 

Then, after integration in momentum and summing in spin, we obtain 


pR 

± = - jjr'te - “) 




eB 

47T 2 


E^- M * + ) 


Z=1 

00 


- 


z=i 


^2 _ ( M +) 2 - (M +) 2 In 

f - (Mr) 2 - (Mr) 2 In 


M + ^ 2 -(M +) 2 
M+ 


m + Vm 2 -(md 2 

Mr 


with 


= V2eBA + m 2 ±77, l ^ 1 


(A16) 


(A17) 


where we have neglected the radiative corrections to the mass for all LL’s, including the LLL. As 
already mentioned, in the strong-field limit, 71 <C 71 for l > 1 (see Fig. 1), which allows us to neglect 
in (IA17I) the contribution of 71 for l > 1. Hence, (1A16I) reduces to 


pR 


eB 

47T 2 


0 (/z-M+) 




_ ra * 1 „ ( 

\ m 

py]y? - (M +) 2 - (M +) 2 hr 

pyjp 2 - (Mr ) 2 - (Mr ) 2 In 


/x±^ 2 -(M+) 2 
M+ 


P+ \/P 2 ~ ( M i ) 2 
"Mf 


- =-* £«(/<-mo 


1—2 


/VM 2 - (Mz ) 2 - (M ,) 2 In 


M ±yV 2 -(Mz ) 2 

Mi 


(A18) 


where 


Mz = \/2eBl + ttz 2 , l > 1 


(A19) 


Expression (IA18I) coincides with Eq. (I5U1) . 
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3. Finite-temperature contribution 

The contribution at finite temperature is given by 



After summing in 77 and <7, and taking into account the corresponding spectra ED, we obtain 

Q%= - dp 3 ln(l + e -0IV'»’! +m2+ ' 1 l)(l + e -MVp 2 3+ m2 -v\) 

p J 0 

dp 3 ln(l + e^ | V / P3+< 2 - +Ml)(i + e -^l V / Pi+< 2 ^l) 
dp 3 ln(l + e _/31 y/pl+ M i 2 +i i \) (l + e-^ 1 ) 

-V / dp 3 ln(l + e - /J| Vp?+ M « 2 +^l)(l + e -/ 3 l\/pi+ Mi2 -^l) ; (A21) 

77 P'S 

Here we are using the notation given in (IA17I) and (IA19I) . 


eB f°° 

2tt 2 /5 J 0 

eB r°° 
2tt 2 /3 Jo 


Appendix B: Thermodynamic potential a la Dittrich 


In this Appendix, we will discuss the steps followed in [23] to calculate the vacuum contribution 
of the thermodynamic potential for charged fermions endowed with AMM in a magnetic field. In 
doing that, we will call attention to some issues we have found in those calculations that explain the 
difference between the strong-field results (1351) and (1551) . 

Let us start by writing Eq.(2.1) of Ref. [23 

(m - ^a^F^ + 7 ■ n)G(a;, y; A) = 5(x - y) (Bl) 

Taking into account that the fermion propagator can be written in term of a scalar propagator, 
A {x,y;A), as 


G(x, y ; A) = (m — - 7 ■ n) A(x, y, A), 


then, we can write Eq. ED in term of A(a;, y; A ), as 

(m - + 7 ■ n)(m - - 7 • n) A (a;, y; A) = 8{x - y), 

which, with the help of the identity [ 7 p ,er M! ,] F A “'n p = 4 i 7 p F pp n p , simplifies to 


(m - 2 - (7 • H) 2 — 2 ^ 7 p E p 


. 11 " 


A {x,y;A) = 6{x - y) 


(B2) 


(B3) 


(B4) 


By comparing Eq. (IB4I) with Eq.(2.3) of Ref. [23[, we can see that the last term inside the squared 
parenthesis in El is missed in Ref. [23| . On the other hand, if we compare the quasiparticle 
spectrum that corresponds to Eq.(2.3) of Ref. 


Po = v\ + 2 leB + [to - T]“ 


(B5) 
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with the very well-known one for charged particles with magnetic moment in a uniform magnetic 
field, we have that the spectrum (IB5I) can be found from (l65l) in the limit 


„ 2 eBl 

1 +— 


(B6) 


which will be only valid at very weak field. Then, the strong-field result (IMl) , which was found after 
the consideration (IB6I) . is not valid. 

To compare our thermodynamical potential with the one reported in E^j . let us rewrite our 
Eq. (1D2I) (see also Eq. (ID4D 1 as 


w _ \ e B\ 

vac 2(2 tt) 2 


+ 


2(2t r) 2 


Applying in Eq. (IB7I) the approximation 

\e.B 


n 


w 


f ~5T y~l ( 2 — Sio)e s ( m +P±+ r ) cos h f 2 sT Jm 2 + p]_ 
Jo s i=o ^ 

l CB l B shlH2sTm) 

we obtain 

/»00 7 

2(2*)* l ^ (Wm) 


(B7) 


1=0 


+ ^pL $ e "' < “’ +T ’ )3il ’ h(2>Tm ) 


Once we perform the summation over all Landau levels, we arrive at 


n w = - 

“ nr- 


\eB\ 

2(2tt)- 


r OO 7 

/ (_£ e -s( m +T ) {coth (eBs) cosh (2sTm) — sinh(2sTm)} 

Jo s 


(B8) 


(B9) 


In [23|, expression (IB9I) was simplified even more considering that T eB and then taking the 
small-7 - limit for the functions cosh and sinh to obtain, after renormalization, the expression (2.10) 

of Ref. El, 


n = - 


2(2t r) 2 


r°° ds_ ( m 2 +r 2 ) 


eBs coth (eBs) — 1 — - (eBs) 2 

O 


(BIO) 


In our calculation, we went beyond the leading term in the expansion at small-7 - of the cosh and sinh 
functions, since there are some extra 7 - -dependent terms that contribute to the divergency. Then, 
in addition to the usual renormalization subtraction carried out in (IB 101) , we had to subtract several 
T-dependent terms (see Eqs. m-m)- Nevertheless, in the leading weak-ficld approximation 
taken up to 0((eB) 4 ) order and in the a 2 correction, the effective potential (1B10I) reduces to the 
one we found in (ID10I) for the weak-field approximation. 


Appendix C: From Ritus to Schwinger propagator 

by using the Schwinger proper 


To perform the summation over all LL’s, let us rewrite Eq. 
time representation 


/»00 _/» 74 

G(x,y) = -i(fl x +m) j ^ ds e is ^- m2+ie) E p (x)H(l)E p (y) 


(Cl) 
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Now, writing the Ritus eigenfunctions explicitly in Eq. m, we obtain 

G(x,y) = -i(fl x + m) [ ds [ ^I!^ e is (Pu- m2 +^) e -i[^o-Vo)po-{x 2 -y 2 )p 2 -(x 3 -y 3 )p 3 ] 

Jo J (2tt) 

OO 

x E e- is(2eBn E N*D n (p x )D n (p y )Il(l)A(a) (C2) 


1=0 


Replacing one parabolic cylinder function by its integral representation in Eq. (IC2I) , and after some 
manipulations, we get 


G{x,y) = -i{fl x +m) J ds J 

/•OO 

xV2eBe* p * / dt 


m +ie) e -i[(xo-yo)po-(x2-y2)p2-(x3-y 3 )p3] 


e -|t 2 girpxt 


“ (-4 r t) n e ~is(2eBl) 

E E 1 -- D n (p y )U(l)A(a) 


r—±1 1 


1=0 cr = ±l 


Performing the sum over a and using the identity 


E^ J D n( z ) = 


,-iz 2 +zt-U 2 


(C3) 


(C4) 


n—0 


on each term of Eq. (IC3I) . we obtain 


G(x,y) = -i(fi x + m)J^ds J ^^ e ^-m 2 + i,) e - 

xV^Be^fl-pl) [ e -^ 2 [l-(e- i( 2 eB 3 ) ) 2 ] e it[Px-p H e- i < 2eB '’)] 

J—OO 


i[(xo-yo)PO-(x2-y2)P2-(x3-V3)P3] 


[A(+) + e~ is2eB A(—) 


(C5) 


where in the last line we have performed the sum over r. 

Once we perform the integration over t and in momentum p , we arrive at the well known result 
for the fermion propagator in an external magnetic field 


G(x,y) = —(yi x +m)$(x,y) - - [ --- e ls ^ m le ^e y h[ 4 tan(es S )(* y A} 

I 67 H J o ssin(ei3s) 

x [e ieSs A(+) + e~ ieBs A(—)] (C 6 ) 

where 

&(x,y) = e iS r lx2 ~ y 2^ Xl +Vl) . (C7) 


is the well known Schwinger phase, which encodes all gauge dependence of the fermion propagator 
(see Ref. j25|). 

It is easy to prove that for the particular gauge A(z) = Bz\z 2 and integrating along a straight-line 
path from y to x the Schwinger phase factor can be written as 

eB f x 

-y(z 2 -J/ 2 )(zi + 2/i) =J dz ■ A(z) (C 8 ) 
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Appendix D: Thermodynamic potential in the weak-field approximation 

1. Vacuum 

By using the proper-time representation, we rewrite Eq. m as 

qW _ l e -^l \ ' I f d -Pll f~is{p^-(m-T) 2 -ie} _ g-isjpj)-(m+<rT) 2 -ie} 

vac ~ 1 a t±Jo 8 J (27r) 3 l 

+ ^(2 _ 5l0 )e- is V 2 \\-^ m ' i+2leB+aT ) ~ u \ L (Dl) 

i=o ) 

where we are using the notation p 2 = pop 0 + P 3 P 3 . 

Note that in (ID 1 1) we added and subtracted the term e ~ ls {p\\~ l ' rn+ljT '> _ze }, and that the factor 
(2 — 5zo) has been conveniently introduced despite once T is present, there is no spin degeneracy for 
higher Landau levels. 

Working in Euclidean space, we now perform in m the integration over p \| and the summation 
over a, obtaining 


q w _ \ e B\ 
vac 2(2 n) 2 


J J e - s ( m2+r2 ) | - sinh( 2 sTm) 

OO 

+ ^(2 - S l0 )e~ 2sleB cosh (2 sTy/m 2 + 2 leB^j 


By using the identity mum 


1=0 


w \ r"t+ ioa ^ 4 _|_x£ 

cosh(cc) = -—: / —e + 4t 

27TZ J-y — ioQ t 2 


with 7 a positive number, we can easily perform the sum over LL’s to get 

\eB\ 


n w = - 

vac 2(27r) 2 


J™ J e - s ( m2+r2 ) | - sinh(2sTro) 


+ 


0r n+ io ° dt t ! s 2 T 2 m? 

27TZ J 'v—ioo t 2 


coth 


S 1 


_ ST2 \ 


t J 


eB 


(D2) 


(D3) 


(D4) 


This particular form of the effective potential allows us to isolate all divergences by making a Taylor 
expansion in powers of ( eBs ) and (T) up to order 0((eB) 4 , T 4 ), it is 




w 


1 


2(2t r) 2 


dS -~ sm2 (\ - sT 2 + \s 2 T 4 ) \eB\s 


x < — 2 smT + 


/»7+ioo 


dt 


27TZ 


7—200 t 2 
3 s 


-e* 1 + 


/ 1 eBs (eBs) 
\eBs + ~~3 45~~ 


s 2 T 2 m 2 s 4 T 4 m 4 
t + 2f 2 

1 1 (eBs ) 2 2(eBs) 4 \ sT 2 eBs 


(eBs) 2 


15 


189 
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f 1 eBs 4(eBs) 3 \ (sT 2 ) 2 (eBs) 
\ (eBs) 3 ~ "l5~ + 189 J f 2 


(D5) 
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Hence, the divergent terms are 


Q w 

^ L div 


ds 


2(2t r) 


— 2 s 2 mTeB 


1 / 7 ? n +io ° dt 


2ni 


' 7 — 200 


t 2 


1 + S (- -11 r 


2 ' s 2 


{eBf 


m 2 T 2 r 4 r 4 r 4N1 


f 


Integrating over t term by term, we arrive at 


nW = _ 

8 tt 2 7 0 . 

where we used the identity 


ds 


1 


1 — 2 ms 2 TeB + sT 2 + 2 s 2 m 2 T 2 — ttS 2 T 4 


(eBs ) 2 


6 


a/ 7T r 4100 dt 


e* = 


2ni 7 7 _ ioo t§+" r(i + n) 



(D7) 


(D 8 ) 


Note that in (ID7I) there are several 7~-de pen dent divergences, similarly to what we also found in the 
strong-held limit. As in the T = 0 case [241 ] . where a vacuum term depending on the renormalized 
mass was subtracted, now we subtract the new magnetized vacuum term, which also depends on the 
particle AMM. At T = 0, the divergent term (1D7I) reduces to that of Ref. [ 13 ]. The last term in the 
RHS of (1D7I) is the usual divergent term proportional to B 2 appearing at T, p = 0 and that can be 
absorbed in the magnetic held renormalization function Z 3 as is usually done in the T = 0 case. 

Then, following the Schwinger procedure of renormalization, we subtract Eq. (lD7l) to Eq. llD2l) . 
to obtain the renormalized effective potential given by 


0) = ~ 2 ( 2 tt ) 2 ,/, $ e s ( m2+r2 )j - |eR|ssinh(2sTm) 


+ 


yJH r +i °° dt t+ £l 


2iri 


^ e t+ '* |eB|scoth |^s ^1 — —) eR 


' 7 — 200 

-2 1 2_2'-r2 


sT 5 


- ( 1 + sT 2 + 2s 2 m 2 T 2 -2ms 2 TeB--s 2 T 4 

6 


(eRs ) 5 


(D9) 


Now e and B have their renormalized values that have absorbed the divergent contribution Cldiv ■ la 
the renormalization procedure we are following, we consider that T is a renormalized parameter, as 
it was considered in the strong-held case in Appendix A, and for the mass in the T = 0 case [241 ]. 
Then, the extracted divergency is considered as part of the renormalization of the vacuum energy 
that now also depends on the AMM. 

In the weak-held expansion, the leading terms of the renormalized potential m , taken up to 
0((eR) 4 ) order and a 2 correction, are 


nZi R \B, 0 , 0 ) 


1 f(e £) 4 

2(2tt) 2 \ 45m 4 


(eR) 2 7^ 1 
3 m 2 J 


(DIO) 


2. Finite density at zero temperature 


The hnite density contribution to the thermodynamical potential in the weak-held approximation 
is given by 


r>w _ \ eB \ 
" ( 2 tt ) 2 


/ dp 3 < (p - £++o )0(p - e++o) + £+ai)0(n-£+ai)> (Dll) 

■’ l cr=±l 1=1 J 
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In the weak-field limit, the sum over all LL’s can be easily done if we rewrite the above equation 
as follows 

^ E J dp 3 ^(p-£++o)0(p-£++o)-{p-£+ao)0(p-£+<To) 

OO 'j 

d" ^ ' (2 ^Zo)(// £j]al)di^p £ 7 ]al') / 


(D12) 


1=0 


Taking into account that the separation between LL’s is given by \f2eB, in the weak-field limit 
the LL’s are so close that we can consider that they form a continuos distribution. Therefore, this 
allows us to use the Euler-Maclaurin formula 32] 


\eB\ 


00 /*00 I o| 

^2(2- S i0 )f(2eBl) ss / dxf{x) + J-y^/(oo) 
z=o ■'° 


OO / T^ \ 

+ £ 7§iT s “ {/‘“-'’(OO) - /‘“-’’(O)} 


/C=l 


(D13) 


where a; is a continuous variable and B 2 k are the Bernoulli numbers with Bq = 1, B\ = — i and 
B 2 = \. 

Thus, (1D12I) up to 0((eB) 2 ) can be expressed as 


oHf = - 


( 2 -) 2 
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d 3 p 


<T =±1 


Of 


(A 4 £+crp±)d(p £+crp±) 
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2 J 

\eB\ 2 

12 


/ #3 


(m — e ++o)f A 4 — £ ++o) — (p — £ +ao)0(p — e+o-o) 
m + <tT 




■Q(p - £+ 0 - 0 ) 


where 


£r/ap ± — 




(D14) 


(D15) 


and we replaced x by i n the first term of Eq. (ID14I) . 

Performing the integration over all momenta and Taylor expanding in powers of 7~, the leading 
contribution is given by 


= ~^~2 ( ~ m?pp F + m 4 In { *' r J ^ ) ]0 (/z-m) 


8 tt 2 \ 3' 
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. o I PPf + m In 
47r z V \ m 


4,„ I PF+ P 
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2 -7T 2 \ m 
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(D16) 


where p F = \J P 2 ~ m 2 . 
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3. Finite-temperature 

The finite temperature thermodynamical potential reads 


flT = - 


\eB\ 


(2tt) 2 /3 


/ ( OO 
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Once we perform the sum over rj we get 
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